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ABSTRACT 

We report the discovery by the CoRoT space mission of a transiting brown dwarf orbiting a F7V star with an orbital period of 3.06 days. CoRoT- 
15b has a radius of 1.12+Jj ^ R Jup , a mass of 63.3±4.1 Mj up , and is thus the second transiting companion lying in the theoretical mass domain of 
brown dwarfs. CoRoT-15b is either very young or inflated compared to standard evolution models, a situation similar to that of M-dwarfs stars 
orbiting close to solar-type stars. Spectroscopic constraints and an analysis of the lightcurve favors a spin period between 2.9 and 3.1 days for the 
central star, compatible with a double-synchronisation of the system. 

Key words, brown dwarfs - planetary systems - low-mass - techniques: photometry - techniques: radial velocities - techniques: spectroscopic 



1. Introduction 



The CoRoT space mission (Baglin et al. 2009), in operation 
since the start of 2007 and extended to the end of 2013, is de- 
signed to find transiting exoplanets. A natural product of this 
mission is that any object with size of Jupiter or lower that 
transits its host star can be detected. This includes stellar and 
sub-stellar companions such as M-dwarfs and brown dwarfs 
(BDs). In the mass-radius diagram of transiting companions or- 
biting solar-type stars, there is up until now only one known 
brown-dwarf, CoRoT-3b (Deleuil, et al. 2008), located in the 
ap in mass between planetary and low-mass star companions 
Determination of the physical properties of such objects are 
fundamental to understand the link between the population of 
planets and low-mass stars and to distinguish the formation and 
evolution processes of the two populations. 

We report in this paper the discovery of a new transiting 
brown-dwarf by CoRoT established and characterized thanks to 
ground-based follow-up observations. CoRoT-15b, with an esti- 
mated radius of 1.12 Rj up and an estimated mass of 63.3 Mj up , 
orbits in 3.06 days an F-type dwarf with solar metallicity. 



* The CoRoT space mission, launched on December 27th 2006, 
has been developed and is operated by CNES, with the contribution 
of Austria, Belgium, Brazil , ESA (RSSD and Science Programme), 
Germany and Spain. Observations made with HARPS spectrograph at 
ESO La Silla Observatory (184.C-0639). 

1 After submission of our manuscript, Johnson et al. (20101 reported 
the discovery of the transiting brown dwarf LHS6343C with a radius of 
0.996 ± 0.026 Rj up and a mass of 70.6 ± 2.7 M Jup . 



2. CoRoT observations 

SRa02 was the seventh field observed with CoRoT in the second 
year after its launch. It corresponds to the third short run and was 
located towards the so-called galactic anti-center direction. This 
run started on 2008 October 1 1 and ended on 2008 November 
12, constituting of a total of 31.7 days of almost-continuous ob- 
servations. 

More than 30 multi-transiting candidates for planets were 
found among the 10265 targets of the SRa02 field. About half of 
them were clearly identified as binaries from light-curve analy- 
sis and around tenth of high priority planet-size candidates were 
selected in this short run including SRa02_El_4106 afterwards 
called CoRoT-15. The various ID of this target, including coor- 
dinates and magnitudes are listed in table [T] 



3. CoRoT light curve analysis 

CoRoT-15, with an estimated V magnitude close to 16, is lo- 
cated at the faint end of the stellar population observed by 
CoRoT . In this magnitude domain, the signal is not sufficient 
to split the photometric aperture in three colors and the ex- 
tracted lightcurve, also called white lightcurve, is monochro- 
matic (Llebaria & Guterman 2006). The sampling rate was at 
512 seconds during the whole run and not oversampled to 32 sec- 
onds since this candidate was not identified with the alarm mode 
(Surace et al. 120081 . Fig. Q] shows the lightcurve of CoRoT-15 
delivered by the N2 data levels pipeline. This lightcurve is quite 
noisy and is affected by several high energy particles impacts 
which result in hot pixels as well as possible stellar variability. 
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Table 1. IDs, coordinates and magnitudes. 



CoRoT window ID 


SRa02.El_4106 


CoRoT ID 


221686194 


USNO-B1 ID 


0961-0097866 


2MASS ID 


06282781+0611105 


RA (J2000) 


06:28:27.82 


Dec (12000) 


+06:11:10.47 




16.85 


B2~ 


16.59 


Rl" 


15.47 


R2 fl 


15.43 


I" 


14.83 


J" 


13.801 + 0.026 


H'' 


13.423 ± 0.037 


K b 


13.389 ± 0.050 



" fromUSNO-Bl catalog 
b from 2M ASS catalog 



1.03F 




0.97 E 

4755 4760 4765 4770 4775 4780 
Time (HJD-2450000.0) 



Fig. 1. Raw light curve of CoRoT-15. 

A first analysis of the lightcurve (LC), based on a trapezoidal 
fit to each individual transit, reveals periodic transits with depth 
of 0.68% and a period of 3.0608+0.0008 days. 

Since the CoRoT lightcurve is relatively noisy, no meaning- 
ful limits to either the visible-light albedo of the companion nor 
to its dayside surface temperature could be established. 

From the raw lightcurve, we tried to estimate the stellar ro- 
tation period. We filtered out variations with timescales shorter 
than 15 data points (~2.13h) to reduce sensitivity to the satel- 
lite orbital effects and other short term variations, and we then 
removed the transits. The computed LS periodogram appears to 
be quite noisy and affected by low-level discontinuities in the 
lightcurve. There is tentative evidence of a possible rotational 
modulation at either 2.9, 3.1 or 6.3 days in the light curve, but 
the data does not enable us to estimate the period more precisely 
or to distinguish between these values. 

4. Ground-based observations 

4. 1 . Ground-based Photometric follow-up 

Ground-based photometry was performed with the aim of refin- 
ing the target ephemeris, to verify that none of its closest con- 
taminant stars correspond to an eclipsing binary and to deter- 
mine the contamination from nearby stars inside CoRoT\ pho- 
tometric aperture mask (Deeg et al. (2009). CoRoT-15 was ob- 



served during a transit event on 2010 January 13 with time-series 
photometry at the IAC80 telescope, from HJD 2455209.480 to 
.676. The resultant lightcurve was not sufficiently precise to 
identify the expected transit on the target, due to photometric 
errors introduced by the presence of thin cirrus. However, the 
absence of large brightness variations in the neighboring stars al- 
lowed us to exclude nearby eclipsing binaries as a source of the 
signals that were observed by CoRoT. The contamination factor 
was derived from a measure of the distance and brightness of 
the nearby stars on a subset of these IAC80 R-filter images ob- 
tained with the best seeing of that night (1.7 arcsec). Ten nearby 
stars were identified, with six of them contaminating the CoRoT 
window aperture with a flux level that amounts to 1.9+0.3% of 
the main target flux. We checked furthermore that none of the 
known nearby stars is bright enough to be contaminating eclips- 
ing binaries. The image of the sky around CoRoT-15 is shown in 
Fig. [2] 



1 arcmin 
M ► 




E 



Fig. 2. The image of the sky around CoRoT-15 (star in the cen- 
tre). Left: R-filter image with a resolution of 1.7" taken with the 
IAC80 telescope. Right: Image taken by CoRoT, at the same 
scale and orientation. The jagged outline in its centre is the pho- 
tometric aperture mask; indicated are also CoRoT's x and y im- 
age coordinates and positions of nearby stars which are in the 
Exo-Dat (Deleuil et al 2009) database. 



4.2. Radial velocity follow-up 

Radial velocity (RV) observations of CoRoT-15 were performed 
with the HARPS spectrograph (Mayor et al. 120031) based on the 
3.6-m ESO telescope (Chile) as part of the ESO large program 
184.C-0639 and with the HIRES spectrograph (Vogt et al. [1995b 
based on the 10-m Keck-1 telescope as part of NASA's key sci- 
ence project to support the CoRoT mission. 

HARPS was used with the observing mode obj_AB, with- 
out simultanenous thorium in order to monitor the Moon back- 
ground light on the second fiber. The exposure time was set to 1 
hour. A set of 9 spectra was recorded between November 24th 
2009 and February 21th 2010. We reduced HARPS data and 
computed RVs with the pipeline based on the cross-correlation 
techniques (Baranne et al. 1 19961 Pepe et al.|2002). The signal-to- 
noise ratio (SNR) per pixel at 550 nm is in the range 3 to 7.8 for 
this faint target. It corresponds to the faint end in magnitude for 
HARPS follow-up observations. Radial velocities were obtained 
by weighted cross-correlation with a numerical G2 mask. 

HIRES observations were performed with the red cross- 
disperser and the L-cell to measure RVs. We used the 0.861" 
wide slit that leads to a resolving power of R * 45,000. The 
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Table 2. Radial velocity measurements of CoRoT-15 obtained 
by HARPS and HIRES. BJD is the Barycentric Julian Date. 



CoKoT 15b 



HARPS l KECK 



BJD 

-2400000 


RV 

[kms- 1 ] 


+ lcr 
[kms- 1 ] 


HARPS 3.6-m ESO 


55159.81312 


9.944 


0.308 


55169.77922 


2.107 


0.315 


55235.60498 


2.849 


0.380 


55236.59903 


8.057 


0.392 


55241.55686 


0.119 


0.338 


55243.55544 


-2.422 


0.222 


55246.56091 


-2.684 


0.434 


55247.55392 


-1.318 


0.411 


55248.54372 


10.235 


0.357 


HIRES 


10-m Keck 


55221.76973 


1.084 


0.291 


55221.81404 


0.890 


0.189 


55221.82490 


0.556 


0.289 


55221.89482 


-0.232 


0.220 


55221.90689 


-0.492 


0.216 


55221.97788 


-1.869 


0.353 


55221.99025 


-1.750 


0.147 


55222.98656 


-3.193 


0.340 


55222.99752 


-3.587 


0.226 


55223.73904 


6.527 


0.113 


55224.82036 


1.593 


0.374 


55224.83152 


1.818 


0.303 


55225.04022 


-1.345 


0.515 



contamination of the HIRES spectra by scattered moon light was 
significant for this faint target, but the 7" tall decker allowed us to 
properly correct for the background light. Three spectra without 
the Iodine cell were obtained on 2010 december 2 and january 9. 
These 3 spectra were co-added to serve as stellar template for the 
RV measurements, and to be used for the determination of stellar 
parameters (see Sect. 14.31) . Over a 4-nightrun from 2010 January 
25 to January 28 we have collected 13 spectra of CoRoT-15 with 
the 12-cell. The average SNR of the spectra with the 12-cell range 
from 13 to 22 (per pixel) in the iodine region from 500-620 nm. 
Differential radial velocities were computed using the Austral 
Doppler code (Endl et al. 120001 1. Nine RV measurements were 
made during a transit event on 2010 January 25 but were not sen- 
sitive enough to detect the signature of the Rossiter-McLaughlin 
effect expected to have, for this high rotating star, an amplitude 
of about 100 ms _1 . 

The HARPS and HIRES radial velocities are given in 
Table [2] The two sets of relative radial velocities were simul- 
taneously fitted with a Keplerian model, with the epoch and pe- 
riod of the transit being fixed at the CoRoT value and with an 
adjusted offset between the two different instruments. No sig- 
nificant eccentricity was found and we decided to set it to zero. 
We found a systematic shift in phase using the CoRoT period 
of P=3.0608 days. It comes from the fact that the quite large 
uncertainty on the CoRoT period (69 seconds) may introduce 
after one year a systematic shift of more than 2 hours. When we 
adjust the period with the RVs fixing the transit epoch as de- 
termined from CoRoT lightcurve Tt=54753.5570 + 0.0028, the 
best solution is obtained for P=3.06039 + 0.00014 days and a 
semi-amplitude K-l. 376 + 0.090 kms . The dispersion of the 
residuals is 0.325 kms -1 and the reduced^ 2 is 0.90. The joint 
analysis of the photometric and RV data, presented in Section|5j 
does not change significantly the results. Figure [3] shows all the 
radial velocity measurement after subtracting the RV offset and 
phase folded to the updated orbital period. 
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Fig. 3. Phase-folded radial velocity measurements of CoRoT-15 
with HARPS (dark circle) and HIRES (open triangle) 



4.3. spectral classification 

Three HIRES spectra of CoRoT-15 were acquired without the 
iodine cell. Each spectrum was set in the barycentric rest frame, 
cleaned from cosmic rays and from the moon reflected light. The 
co-addition of these 3 spectra results in a master spectrum cov- 
ering the wavelength range from 4100 A to 7800 A with a SNR 
per element of resolution in the continuum ranging from 20 at 
5300 A up to 70 at 6820 A. The 9 co-added HARPS spectra 
unfortunately did not permit to reach a better SNR. 

From the analysis of a set of isolated lines, we derived a 
vsin/ of 19 + 2 kms -1 . The spectroscopic analysis was car- 
ried out using the same methodology as for the previous CoRoT 
planets and described in details in Bruntt et al. (2010). However, 
the moderate SNR of the master spectrum of this faint target, 
combined to the marked rotational broadening of the spectral 
lines prevented an accurate measurement of the star's photo- 
spheric parameters. The derived stellar parameters are reported 
in Table E 

Following Santos et al. (2002) methodology, we also esti- 
mated the vsin / and an [Fe/H] index from the HARPS cross- 
correlation average parameters (FWHM and surface). Assuming 
a B - V of 0.5, we estimated the vsin/ of the target to be 
16+1 kms -1 an [Fe/H] index close to zero (solar metallicity) 
in agreement with the spectral analysis. 



5. System parameters 

The time span of the CoRoT light curve is relatively short, and 
the RV data was collected on year later. Jointly analysing the 
two datasets therefore yields significantly improved constraints 
on the period P and time of transit centre T tr . To do this, we used 
a Metropolis-Hastings Markov Chain Monte Carlo (MCMC) al- 
gorithm (see appendix 1 of Tegmark et al. 2004] for a general 
description of MCMC algorithms and Winn et al. l2008l and ref- 
erences therein for a detailed description of their application to 
transits). This has the added advantage of yielding full posterior 
probability distributions for the fitted parameters, ensuring that 
the effects the well known degeneracy between the orbital incli- 
nation i and system scale a/R^ (which leads to highly skewed 
probability distributions for these parameters, as well as for the 
radius ratio R c /R^) are properly taken into account in the final 
uncertainties. 
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The light curve was first pre-processed to remove out-of- 
transit variability as follows. Outliers were identified using an 
iterative non-linear filter (see Aigrain et al. 2009 ), and a straight 
line was fitted to the region around each transit. The effect of 
contamination reported in section l4~T1 was taken into account by 
subtracting a constant amount of flux equal to 1.9% of the mean 
flux from the light curve. Each section of the light curve was 
thus normalised, and a visual check was performed to ensure 
that no residual discontinuities affected the preprocessed light 
curve. The photometric uncertainties were then estimated from 
the out-of-transit scatter in the preprocessed light curve section 
around each transit. The light curve was modeled using the for- 
malism of Mandel & Agol d2002l >. Given the relatively low SNR 
of the transits, we opted to fix the quadratic limb-darkening pa- 
rameters u a and Uh at the values given by Sing et al. (2010J for 
the star's effective temperature, gravity and metallicity (adopted 
values: u a = 0.32, u/, = 0.30). The RV data were modelled us- 
ing a Keplerian orbit with eccentricity fixed at zero, since the 
data show no evidence for a significant eccentricity. The rela- 
tive zero-point of the HARPS and HIRES velocities, 6Vq, was 
allowed to vary freely. The parameters of the MCMC were thus 
P, T ti , R c /R+, a/R*, the radial velocity semi amplitude K, the 
systemic radial velocity V sys and SVq. 

After an initial chain of 10 5 steps to adjust the MCMC step 
sizes for each parameter, we ran 10 MCMC chains of 10 5 steps, 
each with different starting points. The convergence of the chains 
was checked using the Geldman-Rubin statistic (Geldman & 
Rubin [19921 Brooks & Geldman 1997). The chains were then 
combined (after discarding the first 10% of each chain, where 
the MCMC is settling from its starting point) to produce poste- 
rior probability distributions for each parameters. We report in 
Table Ejthe median of the probability distribution for each pa- 
rameter To estimate uncertainties for each parameter, we com- 
puted the range of values which encloses 68.5% of the proba- 
bility distribution (rejecting 16.25% at each extremum). Our un- 
certainties thus correspond to 68.5% confidence intervals, just 
as classical 1-cr uncertainties do for a Gaussian distribution. The 
best-fit transit model is shown superimposed on the folded light 
curve in Figure|4] To highlight the correlations between b, R p /R* 
and a/R+, and explain the rather large resulting uncertainties, we 
also show the posterior probability distributions and 2-D pro- 
jections of the combined MCMC chain for these parameters in 
Figure 

We used the photospheric parameters from spectral analysis 
and the stellar density derived from the transit modeling to de- 
termine the star's fundamental parameters in the (T e ff, /R+) 
space. Using STAREVOL evolution tracks (Palacios, private 
com.), we find the stellar mass to be M* = 1.32 + 0.10 Mq and 
the stellar radius R* = 1.46 Rq, with an age in the range 
1.14-3.35 Gyr. This infers a surface gravity of log g - 4.23 *02Q> 
in good agreement with the spectroscopic value. 

Calculations using CESA M (Morel & Lebreton 20081 see 
also Guillot & Havel 1201 Oi l confirm these solutions. The age 
constraints, 1 .9+1 .7 Gyr, are however extremely weak, and yield 




Fig. 4. Folded, detrended light curve of CoRoT-15 (top), show- 
ing the best fit transit model (red solid line) and residuals (bot- 
tom). 



0.0 0.3 0.6 0.07 0.08 0.09 




0.0 0.3 0.6 0.07 0.08 0.09 

b R p /R t 



Fig. 5. Selected posterior probability distributions and two- 
dimensional correlations for the transit fit. The panels along the 
diagonal show single-parameter posterior probability distribu- 
tions for b, R p /R+ and a/R^. The red, blue, and dashed blue 
vertical lines indicate the location of the best-fit, median and 
the limits of the 68.5% confidence interval for each parameter. 
The off-diagonal panels show, for each pair of parameters, scat- 
ter plots of 1000 points randomly selected from the combined 
MCMC chain. The density of the points approximates the joint 
posterior probability distribution. 



2 The choice of which statistic to report is a somewhat tricky one. 
When the distributions are (close to) Gaussian, the median, most prob- 
able and best-fit values coincide. When the distributions are skewed, 
as in the case of b for example, the median, most probable and best- 
fit values can differ signficiantly. Whilst the best fit value maximises 
the merit function for the particular dataset being analysed, it has no 
special physical meaning. We adopt the value which divides the proba- 
bility distribution in half, namely the median, as it is arguably the most 
physically meaningful. 



possible pre-main sequence solutions with extremely young 
ages. 

We derived for the transiting companion M c =63. 3+4.1 Mj up 
andR e =1.12^°R Jup . 
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Table 3. Star and companion parameters. 



Parameters 



Value 



Transit epoch T„ [HJD] 
Orbital period P [days] 
Transit duration d„ [h] 

Orbital eccentricity e 
Semi-amplitude A* [ km s -1 ] 
Systemic velocity V§_h arps [kms" 1 ] 
Systemic velocity Voj,, reJ [kms -1 ] 

Radius ratio R c /R+ 

Scaled semi-major axis a//?* 

Impact parameter b 

Ml l3 /R+ [solar units] 
Stellar density p* [g cm -3 ] 
Inclination i [deg] 

Effective temperature T e ff\K] 
Surface gravity log g [dex] 
Metallicity [Fe/H] [dex] 
Rotationl velocity v sin i [ km s~'] 
Spectral type 

Star mass [M Q ] 
Star radius [R Q ] 
Distance of the system [pc] 

Orbital semi-major axis a [AU] 
Companion mass M c [Mj up ] 
Companion radius /? t .[Rj U p] 
Companion density p c [g era -3 ] 
Equilibrium temperature T p e " [K] 



2454753.5608±0.0011 

3.06036±0.00003 

3.24±0.1 

(fixed) 
7.36+0.11 
2.23+0.11 
1.09+0.07 



n msfi+o.0039 

U - U/8S -0.0029 



6.68 + "- 49 ' 
0.38+* 



75 +0 05 
60+ 8 ' 3 
86.7111 



6350±200 
4.3±0.2 
0.1 ±0.2 
19±2 
F7V 

1.32±0.12 
1 46 +031 
1270±300 



045 +0 014 
63.3±4.1 

1 1 9+0.30 

,=0.15 

59 +37 

J '-32 

fl20 



1740 4 



!■)() 



6. Discussion and Conclusion 

CoRoT-15b is one of the rare transiting companion that lies in 
the theoretical mass domain of brown dwarfs (13-75 MJup, if 
one adopts the present IAU convention). Contrary to CoRoT-3b 
(Deleuil et al. 2008 ) that is located in the overlapping region be- 
tween the massive planet and the brown-dwarf domain, CoRoT- 
15b is well in the mass domain of BDs. Expanding a bit the 
mass domain, one can easily include in this ensemble the high 
mass "planets" (M > 10 M Ju p ) XO -3b (Johns-Krull et al.,2008) 
and WASP- 18b (Hellier et al. 2009}, and in the M-dwarf regime, 
OGLE-TR-122b (Pont et al. I2005al) -123b (Pont et ai. l2(f)6t - 
106b (Pont et al. 12005bK and HAT-TR-205-013 (Beatty et al. 
120071 ). Interestingly, all these objects are found to orbit F-type 
stars (see also Deleuil et al. 2008), with one exception: OGLE- 
TR-122b orbits a G-type dwarf but has a much longer orbital 
period (7.3 days compared to less than 4.3 days for all other ob- 
jects). 

Early- and mid-F-type dwarfs have the particularity of be- 
ing fast rotators, independently of their age (Nordstrom et al. 
119971) . a consequence of a small or inexistent outer convective 
zone, weak stellar winds, and reduced losses of angular momen- 
tum. The tides raised on a star by its close-in companion (planet, 
brown dwarf or star) have long been known to pause a threat to 
its survival (e.g. Patzold & Rauer 2002). This is true when the 
star's spin is slower than the orbital period of the companion, 
a common situation for close-in exoplanets. However, massive- 
enough companions have the possibility of spinning-up the star 
and may escape engulfment if the total angular momentum of 
the system is above a critical value (Levrard et al. 2009). Even 



in that case however, magnetic braking in the central star (e.g. 
see Barker & Ogilvie 2009) will lead to a loss of angular mo- 
mentum that will be transferred to the orbit of the companion 
through tides and lead to orbital decay. We thus propose that 
close-in massive planets, brown dwarf or M-dwarf can survive 
when orbiting early or mid F-type dwarfs but that they tend to 
be engulfed by G-type (or late F-type) dwarfs. In the case of 
CoRoT-15, we thus expect that the star should be above ~ 1.25 
M Q to avoid efficient spin-down, and that the system should be 
at or close to double-synchronisation (i.e. the spin period of the 
star should be close to the orbital period of its companion). 
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Fig. 6. Masses and radii of eclipsing brown-dwarfs and low- 
mass stars (circles with error bars, as labeled) compared to the- 
oretical mass-radius relations (lines). The lines correspond to 
isochrones of 0.5 (upper orange lines) and 5 Ga (lower blue 
lines), respectively. The dashed lines are calculated for isolated 
brown-dwarfs and low-mass stars. The plain lines include the 
effect of irradiation with T eq — 1800K. The dash-dotted lines 
include irradiation and account for a 50% coverage of the photo- 
sphere with low-temperature spots (see text). A 5 Ga isochrone 
for isolated brown dwarfs/M stars from Baraffe et al. (2003) is 
shown for comparison (dotted line). 



It is interesting to see that given the vsin/ and stellar ra- 
dius determinations, the projected spin period of the central star 
is P/sini* - 3.9^'j days. An LS periodogram shows the pres- 
ence of many peaks possibly due to low-level discontinuities in 
the lightcurve. The most robust peak compatible with the vsin/ 
determination lies between 0.32 and 0.34 cycles/day, and may 
thus be linked to a stellar spin period between 2.9 and 3.1 days. 
The CoRoT-15 system thus appears to be indeed close to double- 
synchronous. Further observations of the system and in particu- 
lar a precise determination of the stellar spin period would be 
a powerful mean of understanding the dynamical evolution of 
this system. Coupled to studies of similar systems, this may also 
yield strong constraints on the on the tidal dissipation factor in 
F-type dwarfs. 

CoRoT-15b is also extremely interesting for its size in com- 
parison with other objects in this mass range, and of evolution 
tracks for hydrogen-helium brown dwarfs and stars. Figure [6] 
shows that it appears inflated compared to standard evolution 
tracks for these kind of objects (Baraffe et al. 2003), although 
it may be compatible with a young age if the true size is at the 
lower-end of the one inferred from our measurements. However, 
we notice that the same problem arises for OGLE-TR-123b, 
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OGLE-TR-106b and, but to a lesser extent, HAT-TR-205-013.0 
The two other known brown dwarfs with direct radius measure- 
ments, discovered in the 2MASS J05352184-0546085 eclipsing 
binary system (Stassun et al. 2006), have very large radii (5.0 
and 6.5 Rj up ) but related to the very young age of the system (~ 
1 Myr), still in the earliest stages of gravitational contraction. 

In order to examine possible solutions to this puzzle (other 
than a systematic overestimation of the inferred sizes for the 
systems known thus far), we calculate evolution tracks using 
CEPAM (Guillot & Morel 1 19951) . but adding the dominant ther- 
monuclear reaction cycle, namely the pp-chain (see Burrows & 
Liebert fT993l l. The atmospheric boundary condition is adjusted 
to the Baraffe et al. (2003 ) evolution tracks, using the analytical 
solution of Guillot (2010), and values of the thermal and visible 
mean opacities, k± - 0.04 cm 2 g~' and kv = 0.024 cm 2 g _1 . The 
model shows that irradiation effects, although quite significant 
for Jupiter-like planets have rather small consequences in the 
brown dwarf regime, and become completely negligible in the 
stellar regime. We also test the possibility that these inflated sizes 
may be explained by the presence of cold spots on the brown 
dwarf, similarly to a mechanism proposed to explain that M- 
type star in close-in binaries also appear inflated (Chabrier et al. 
2007 ). As shown in Fig. [(J this mechanism works only in com- 
bination with a young age for the system, and with a rather large 
fraction (~50%) of the photosphere covered with spots. An alter- 
native possibility is that irradiated atmospheres are much more 
opaque than usually thought, possibly a consequence of photo- 
chemistry and disequilibrium chemistry. A solar metallicity was 
assumed for the brown-dwarf models displayed in Fig. 6. We 
tested the effect of metallicity on the mean molecular mass and 
the opacity but did not find significant change in the radius. 

In any case, this shows that CoRoT-15b is a crucial object 
to understand both the dynamical and physical evolution of gi- 
ant planets, brown-dwarfs and low-mass stars. Further observa- 
tions aiming at obtaining more accurate spectra of the star would 
be highly desirable. Although this is a challenging measurement 
given the faintness of the target, measurement of secondary tran- 
sits in the infrared would be extremely interesting because they 
would inform us on this rare heavily irradiated brown dwarf at- 
mosphere. 
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